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I. REVIEW OF OUR PREVIOUS WORK 

In jl| , the factorization of the Hamiltonian of the quantum harmonic oscillator has been 
effected by means of the pair of non-mutually adjoint operators 



B- 



-L(a.-W-| + ^W). (la) 

We notice that these operators coincide with the standard creation and annihilation oper- 
ators a and a* when a{x) = 1 and (3{x) = x and with the operators introduced by Mielnik 
{2! b and b* when a{x) = 1 and /3{x) = x -\ — ^p^^ ^ where F{x) = e~^^dt is related to the 
error function E(x) through F(x) = ^E(x) while 71 is the integration constant parameter 
occurring through the integration of the Riccati equation. The factorization = H +\ 

corresponding to the eigenvalue problem Hip = Xip of the quantum harmonic oscillator is 
effectively written as 
1 



B~B^ 

2 



ap — H hp 



dx^ \ a J dx a 



so that a{x) and /3{x) should fulfill the following coupled equations 

a' + Pa"^ - /3 = 0, (3a) 

(3' + a/32 = (1 + x^) a. (3b) 

By decoupling these equations one gets the Riccati equation: 

d f . f 1,2 /A^ 

= 1 + X . (4) 



dx \ a J \a 

In a particular solution has been obtained by taking /3(x)/a(x) = x which leads to the 
following simple form of the functions a and /3 

1 X 

as{x) = , =, Psix) = , =, S = constant. (5) 

V 1 + (5e-^ VI + 5e^^ 

To avoid singularities, the parameter 6 should be in the range — 1 < 5 < 00. 

The goal of this Letter is to obtain a complete general solution of the above Riccati 
equation (jlj) which occurs in the alternative factorization based on operators which are not 
mutually adjoint. This solution enables one, as shown in section II, to construct a two- 
parameter family of self adjoint operators from which the standard harmonic oscillator, the 



one-parameter oscillators introduced by Mielnik, and the Hermite operator can be obtained 
in particular cases. In section III, we discuss another possibility of factorization of the 
quantum harmonic oscillator, apart from the one addressed in the previous section. Finally, 
we summarize our work in section IV. 



II. TWO-PARAMETER SELF-ADJOINT OPERATOR FROM ALTERNATIVE 
FACTORIZATION 

Here, we are interested in the most general solution of @. Since this is just the Riccati 
equation for the standard harmonic oscillator, the solution is well known 

— = X ^ . (6) 

a 71 + F(x) ^ ^ 

Using (3{x) from ([6]) in fl3a|) we get 

a' + a (a"^ - l) ( x H ^ , 1=0, (7) 

^ ^ \ 71 + F(x) y ' ^ ^ 

which is a Bernoulli equation that can be solved by the method of separation of variables. 
The general solution reads: 

yi -72(71 + F(x))-'e-^ 

where 72 is a Bernoulli integration constant. In the following, we choose the solution with 
the plus sign because due to eq. ([6]) the other sign cancels out in operators product ([2]). 
Thus: 

'^7172 ('^) / ) (9) 



y^(7i + F(x))^-72e- 



and 

/57.72(x) = ; + + _ ^^^^ 

yi -72(71 + F(x)) 'e-' 
The subindices could be occasionally omitted henceforth. Since the operators B~ and 
define the Hamiltonian, they should be well behaved for any x. This depends on the value 
of the parameters 71 and 72 in a^j^^^^x) and /37i72(^) discussed further. It is also worth 
noting that the product (^^) in B~B^ looks similar to the Hartle-Hawking factor 

ordering js] in SUSY quantum cosmology {4]. 



From the expressions of the operators ([T]) and of (Q and (fTOl) . it is clear that we have to 
determine the values of the parameters in such a way that not either zero or not 

defined. From the numerator of and taking into account that < F{x) < \/7r/2, 

the condition |7i| > ^ does not allow a^^^2(^) to become zero. This is exactly Mielnik's 
condition on his parameter 7. In addition, if 72 < the denominator cannot be zero. 
On the other hand, to see what happens in the case 72 > is not as simple because the 
resulting inequality is transcendental. The fact that (71 + F(x))^ — 726"^^ > leads to: 
(71 + F(x))^ > 726"^^. Since 726"^^ < 72 for all x we will consider only those values of 71 
for which (71 + F(x))^ > 72. Taking the maximum and minimum allowed values of the error 
function in this inequality we obtain ^71 + > 72 and ^71 — ^ j > 72, respectively. 

Their product reads |7i — f | > 72, where we have eliminated the modulus in the right 
hand side because we consider only the case 72 > 0. However, a numerical approach shows 
that an inequality which better forbids the appearance of singularities in the transcendental 
equation is 72 < 7^ — 1, and therefore, the two parameters must satisfy the inequalities 
I71I > y/n/2, and 72 < 7? - 1 or 72 < 0. 

The various factorizations can be now obtained from the general factorization as follows: 
The functions as{x) and (3s{x) found in l| are obtained if we impose the conditions: ^ = —5 
and I71I 00. For the case of the factorization introduced by Mielnik it is enough to take 
72 = and 7i = 7. Finally, the case of the standard factorization is obtained if we have 
•^7172 (^) = 1 /3^^^2(x) = X, such that 72 = and |7i| tends to infinity. 

In Figure dl we can see the geometric positions in the parameter space defining the various 
factorizations, with the x-axis defined by 71 = I/71 and the y-axis defined by 72. The reason 
to use the inverse of the parameter introduced by Mielnik is to provide a more convenient 
display of the factorizations in the two-parameter space, but perhaps it is interesting to note 
that it has been used before in other SUSY developments [s], 0]. Thus, in fig. [T]one gets: (i) 
The standard factorization 7f ^ — )■ ±00 or 71 — )■ 0, 72 = corresponds to the origin in this 
new parameter space; (ii) The factorization introduced by Mielnik: 72 = 0, on the horizontal 
axis between (— ■^); (iii) The factorizations with the help of operators which are not 
mutually adjoint in |l| are encountered when one moves along the vertical axis. The case 
5 — — 1 is obtained when 72 — )■ +00. If 5 > 0, we move along the negative vertical axis with 
the Hermite equations corresponding to the limit 72 — ?■ —00. 

We also display several graphics of a^^^^{x) with the goal to show how it chang 
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function of the parameters 71 and 72 in Fig. |2l In the case of the curves that he above the 



hne a 



71 72 



1, one can see that the peaks pass from the first quadrant to the second one or 



viceversa when the sign of 71 is changed, see, e.g., the curves corresponding to the pairs of 
parameters (1.5, 1.249) and (-1.5, 1.249). If we take values very close to the dashed (blue) 
curve in Fig. [H then the peak of a^^y^ [x) becomes bigger at higher values of the parameters. 
All these curves have been scaled in order to have all of them in a single plot. Notice that 



for values close to 71 = ±-\/7r/2, for 72 < 0, below the line a. 



^7172 



1. On the 



other hand, some plots of fi.^^^^[x) for the same values of the parameters are shown in Fig. 
El Since asymptotically a^^^^{x) — )■ 1 and the second term in the right hand side of eq. ([6]) 
goes to zero, we can see the asymptotic bisectrix behaviour of this coefficient. 

We now deal with the reversed product operator B^B~ . It is more convenient to add \ 



to relate it to the Dirac factorization of the harmonic oscillator. Then, L 



71,2 



leads to: 



L71.2 - 2 



272a/3e" 



d 



dx^ (^^ + F(x))^ dx 
Introducing the functions: 



+ (l + a^) (3^ - (l + x^) + 1 



fill 



(12) 



where il^n{x) are the eigenfunctions of the harmonic oscillator and applying the L^^ ^ operator 



I w^i^^ 

^71,2 -"^n+l 



B+B- + - ) = B' 



B-B+ + - ) ^„ = (A„ 



1) , (13) 



where A„ = n + | are the harmonic oscillator eigenvalues, one can see that they are the 
eigenfunctions of the L^^ ^ operator, but without including the Hq^'^ function. However, the 



latter one can be introduced in the usual SUSY manner asking that L^^ 

B+B 



71,2 



AoH^^", or 



- 1 H^^'^ 



2 ° 



This leads to: 



which requires 



B-K 



71,2 



1 

7! 



a. 



7172 



B+B-Yi^'^ 



1 d_ 



0, 



X + 



71 



Yix) 



(14) 



(15) 



(16) 
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The solution of this equation is: 



(71 + F(x))^ - 726- 



(17) 



where V'oAf(a^) 



■00 



is the modulated zero mode obtained by Mielnik. In fig. H] 



l7i + F(x) 

various plots of Hq^'^(x) are displayed for some values of the parameters. When 72 — +00, 
these functions become more and more singular. 

Since in general any second-order differential equation Pu" + Qu' + Ru + Am = can be 
transformed to the self-adjoint form ^ {j)u') + qu + \u{x)u = 0, multiplying it by the factor 
(1/P)exp (^J^ {Q/ P)dx) , we find in our case that this factor is —2a~^{x) which leads to the 
following eigenvalue equation: 



L,,,,Hr (x) + A„a;,,,(x)Hr (x) = 0, 



(18) 



where 



726" 



;7i + F(x 




726" 



:7i + F(a;)) 



2 ' 



+ x^ 



I3m = x + 



1 - 



726 



:7i + F(x) 



Pi 



F'(x) 



71 + F(x) 



(19) 



is a new self-adjoint harmonic oscillator operator with cj. 



71,2 



X 



^•^7172 (^) weight func- 



tion, which according to the general Sturm-Liouville theory should be strictly positive except 



possibly at isolate points where u. 



71,2 







7|. 



A convenient expression for the eigenfunctions of the operator L^^ ^ can be obtained by 
writing in terms of a and a* because applied to the functions tpn convert them into 
the eigenfunctions B.'J^^-'^ (x) of 3: 



^^TiTs I ^ I 



71 + F(a;) 

This implies the following relationship between the two sets of eigenfunctions: 



(20) 



r7i.2 , 



^7172 ^ 



Vn + 1 ipn+i{x) + 



i'nix) 



(21) 



71 + F(x) 

These eigenfunctions are orthogonal because their construction is performed according to the 
Sturm-Liouville theory and correspond to the same equidistant harmonic oscillator spectrum 
because is isospectral to the harmonic oscillator Hamiltonian ( IT8l) . 
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Figures H] and [5] display the plots of the ground-state eigenfunction Hq^'^ (x) and first 
excited state eigenfunction H]^^'^(x) for several representative cases of the parameters 71 and 
72. From the figure [Hone can infer that 71 acts as a shift parameter for the maximum of 
the wavefunction along the x-axis. The physical interpretation of 72 is not as easy and clear 
but it has to do with the distortion of the shape of the eigenfunctions. 

It is now possible to comment on why in Fig. [H we indicate a region for the Hermite 
equation. The reason is the following. Assume that |7i| — t- 00 and at the same time the 
quotient —72/71 = 5 is big but still a finite number, then the self-adjoint operator ( IT9|) takes 
the form: 

and moreover 

A„w^, ,(x) ^ + 25e-^' = (2n + 1) ^e"^'. (23) 

On the other hand, 

H;^^-^(x)^v^r^H„(x) , (24) 

where B.n{x) are the Hermite polynomials. Therefore, in this approximation, eq. (|T8|) turns 
into 

~ ^"^^ + ^''^") = °' ^^^^ 
which is Hermite's equation. The way in which the functions H^^'^(x) for n = and n = 1 
tend to the Hermite polynomials can be seen in figures. H] and [51 

Similarly to the Mielnik case, the functions HJ^^'^(x) do not admit first-order creation 
operators. To show this fact suppose that the creation and annihilation operators for H^^'^ (x) 
are and A, respectively. Then one has: 

A^Rl-^^ix) = c,R:i\ix), A Rl^^^ix) = c.iQlMx), (26) 

where Ci and C2 are constants. Using fll2p one can realize that applying the operator we 
get: 

B^Rl-'ix) = B-B+ij^^i{x) = (^i/ + ^„„i(x) = n V'„-i(x). (27) 
Applying now the standard creation operator a*, leads to: 

a*B-}il''^{x) = na>„_i(x) = ny/Eiprrix). (28) 
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Finally, by applying we go from ipnix) to }iJ^J^i{x): 

B+a*B-Iil''^{x) = ny/^B+tPn{x) = ny^El'^\{x). (29) 

Therefore the creation operator is defined as: 

A+ = B+a*B- (30) 

and one also gets Ci = n^/n. Since B~^, a* and B~ are linear first-order operators A'^ will 
be a differential operator of third order. 

To calculate the annihilation operator A, we use the results from ( !27l) to which we apply 
the annihilation operator a: 

ai?~HJ^^'^(x) = naipn~i{x) = n\/n — 1 ?/^„_2(x). (31) 

Finally, applying the operator 5+, one gets 

B^aB-Kli^\x) = nVn^B+tljn-2{x) = nV^T^ H^'4(x), (32) 

wherefrom we get the form of the A operator as, 

A = B+aB-. (33) 

Similarly to the operator A~^, the operator A is a differential operator of third order and 
C2 = n\/n — 1. 

Similarly to the factorization introduced by Mielnik, the function of subindex zero cannot 
be obtained using the ^4+ and A pair only. In this case, they transform the functions H^^'^ {x) 
into the functions Hj^!^\(x) or HJ^^'^(x) into B.'J^^^{x) but the constants Ci and C2 are different 
of those obtained when the operators a and a* are applied to ipn- 



III. CAN ONE HAVE OTHER FACTORIZATIONS ? 

One of the obvious questions in the case of the factorization introduced by Mielnik is 
what happens if instead of requiring bb* = aa* one imposes b*b = a*a which looks natural 
because aa* + - = a*a — - = H. However, factorizing in this manner leads to 0(x) = 

^7 — e*^ dt^ e^^ , which is a monotonically increasing function which is singular at the 
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particular value of x where the integral equals 7. Does the same problem occur for our 
factorization ? 

Suppose now that instead of requiring B^B^ = if + | as defined in ([T]) we ask that 

B+B- =H -\. 

We again can get a Riccati equation similar to (jlj) but this time for the product a/3. It 
is then easy to obtain the factorizing coefficients in the following explicit form 



a^,^,{x) = ±Jl + ^ . = ± (34) 

{^i - e* dt) /(^^ _ J- ^t^^^Y + ^^gx^ 



where ki and k,2 are constants. However, one can see that there are problems with the 
continuity of the functions, namely ki — e^^dt could be zero for special values of x. For 
this reason, we present only the particular case 



a{x)f3{x) = X . (35) 



Thus, 



0^,3 (x) = V 1 + 736^' , = . ^ ^ . (36) 

To have B~^ and B~ well defined along the full x axis we require 73 > 0. We notice that in this 
case there is no particular value for the constant that can lead to the factorization introduced 
by Mielnik or to the 5-parameter non- mutually-adjoint factorization. This happens because 
we have started with the different condition B^B^ = H — ^ and we use a particular solution 
of a{x)P{x). However, we can go to the standard factorization by making 73 = 0. 
We now define the operator L^g = B^ B^ — i, or explicitly 

„2 



~ _ _1 _ 733:e^ d_ 1 + x^ + 736^' _ 1 

2 dx^ 1 + 736-' rfx 2(1+736-')' 2 ^ 



and its eigenf unctions H^^ (x) as follows 

Rf{x) = B-4j^+,{x) . (38) 

Thus 

l,3h;^^ = (^B-B+ - (B-V'n+i) = B- {b+b- - = XnK'- m 

Differently from the factorization used in the previous section, in this case it is not necessary 
to calculate Hq^ because it is defined in (|38|) . Moreover B^ = a~^{x){d/dx + x) = a~^{x)a, 
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and therefore the unnormahzed eigenfunctions should be H^^ = a~^{x)ipn{x)- Following the 
analysis performed in the previous section, we will seek the self-adjoint form of the second- 
order differential operator (157|) . It is easy to see that we have to multiply by the following 
factor 




X 



«2 



-2exp(2/ -^^dx] =-2«2(x) = -2(l + 73e^'). (40) 

Under multiplication, ( P7|) shows that the functions H^^ are orthogonal by construction 
leading to the following eigenvalue problem 

l^3h;^3 + ^^^^^(^)h;^3^o, (41) 

where 

= V + ) ir^ + T. + 1 + ;e-' — 

is the one-parameter self-adjoint operator with the weight function O0y^{x) = 2 ^1 + 736^^^ 
which is isospectral to the quantum harmonic oscillator operator, obtained in the limit 
73 0. Interestingly, in the large limit 73 > 1, one can obtain the following differential 
equation 



d ^ d ^ , 

— + 2x— + 2(n + r 

ctx^ ax 



Gn{x) = , (43) 



which differs from the Hermite equation only in the sign in front of the first derivative. 
The corresponding eigenfunctions are of the quantum oscillator type, but with the Gaussian 
factor of double width, i.e., Gn{x) = Cne~^^ Hn{x). 

IV. CONCLUSION 

This work contains the generalization of the factorization procedure introduced by Mielnik 



for the quantum harmonic oscillator performed with a pair o 



non-mutually adjoint factoriza- 



and leads to a two-parameter 



tion operators. It is also an extension of the previous paper 
self-adjoint second-order operator which contains the standard harmonic oscillators, the one- 
parameter harmonic oscillators introduced by Mielnik, and even the Hermite operator for 
special values of its parameters. Both parameters are essentially integration constants, one 
coming from the Riccati equation and the other from the Bernoulli equation which occur 
in the factorization procedure. The Riccati parameter is a shift parameter controlling the 
position of the maxima of the eigenfunctions along the x-axis, while the Bernoulli parameter 
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is related to the shape of the eigenf unctions. In the final part, we have also introduced a 
Bernoulli parameter case which is different of the one discussed in jl| . Moreover, these self- 
adjoint operators are of the effective mass Schrodinger type with position-dependent mass 
which are known to have important applications 8|. 
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FIG. 1: Plot of the allowed region in the parameter space; in particular, we emphasize the values of 



the parameters that lead to t 
the factorization obtained in 



le factorizations studied in this work: the vertical axis corresponds to 
1 

whereas the horizontal axis corresponds to Mielnik's factorization. 
The Hermite equation and the delta-like singular Hq^'^ are found in the extremes of the vertical 
axis. 
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10 



8 - 




FIG. 2: Plots of the functions a{x) for different values of the parameters (71, 72)- The highest peak 
in this panel, corresponding to the parameter pair (71,71 — 10) for 71 = 10^, is also of the smallest 
width. 
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(1.5, 1.249) — 


— (10, 96) — 


(-1, -0.6) 


(10^999990) 


(100, -100) 


(0.9, -10^) 


(-1.5, 1.249) 







FIG. 3: Plots of I3{x) for the same values of the parameters as used in Fig. [2j The difference from 
the standard /3(x) = x plot shrinks at high values of the parameters. 



14 



2 



1.5 - 




(0.9, -10" 
(10^ 999990) 



(100, -10^^) 
(10, 96) 



(10^0) 



(100, -100) 
(-0.9, -10"^) 



FIG. 4: The function Hq^'^(x) for different values of the pair of parameters. Notice the symmetry 
of the red and cyan curves, which differ only in the sign of 71. For all non-zero 72, these functions 
are not Gaussian. 
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FIG. 5: The functions H^^'^(x) for different values of the parameter pairs. 
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(-0.9, -10 ^) 
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FIG. 6: Same as in the previous figure but for tlie functions H2^'^(x). Figure that will not be 
included in the published version. 
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(0.9, -10"^) 

(lolo) 



(100, -10^") (10^ 999990) 



FIG. 7: Same as in the previous figures but for tlie functions H3^'^(x). Figure that wiU not be 
included in the pubhshed version. 
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